Small size boundary effects on two-pion interferometry 
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Bose- Einstein correlations of two identically charged pions 
are derived when these particles, the most abundantly pro- 
duced in relativistic heavy ion collisions, are confined in finite 
volumes. Boundary effects on single-pion spectrum are also 
studied. Numerical results emphasize that conventional for- 
mulation usually adopted to describe two-pion interferometry 
should not be used when the source size is small, since this is 
the most sensitive case to boundary effects. Specific examples 
are considered for better illustration. 



I. INTRODUCTION 

It is generally expected that high energy heavy-ion col- 
lisions may provide the tools to probe the existence of 
a new phase of matter of strongly interaction particles, 
the quark-gluon plasma (QGP), at high temperature and 
high baryon density |jj . The hope of discovering the QGP 
in high energy heavy-ion collisions is to some extent con- 
nected to the possibility of measuring the geometrical 
sizes of the emission region of secondary particles. An im- 
portant tool for accomplishing such size measurements is 
the so-called Hanbury-Brown-Twiss (HBT) interferom- 
etry M. This method was originally proposed in the 
50's for measuring stellar radii but, shortly afterwards, 
it was discovered j| that a similar procedure could also 
be applied to high energy collisions for determining the 
dimensions of pion emitting sources. This method has 
extensively been developed, improved, and better under- 
stood since the pioneering times Gj . 

Differently from the stellar case, however, where the 
dimensions are indeed immense, in the subatomic level 
the effects associated to the small sizes of the particle 
emitters and their boundaries may have an important 
role. Indeed, already in the well-known paper by Gyu- 
lassy, Kauffmann, and Wilson M, and more recently, in 
Rcf . p-pdJI , effects of source finiteness on particle spectra 
and correlation functions were considered, although the 
conclusions of some of them were somewhat contradic- 
tory. For example, the low transverse momentum region 
of Ref. [plpTj is shown to be enhanced with respect to the 
Bose-Einstein distribution. However, this enhancement 
was not observed in other references quoted above. As we 
shall see later, in agreement with results of Ref. |7],p|,p0[ , 
a depletion in the low momentum region is observed in- 
stead. This apparent discrepancy may be explained by 



both the form chosen for the density matrix in Ref. ^"Tj , 
and by the full field theoretical approach adopted there. 
However, the inherent difficulties of that approach are 
enormous, and the simpler treatment discussed in the 
present paper already sheds light to the relevant points 
of the problem. 



The approach suggested in Ref. MM seemed appeal- 
ing for the following reasons. First, it considers that in 
an ultra-relativistic nucleus- nucleus collisions pions are 
the most abundant produced particles, being emitted 
at freeze-out temperatures around 0.1-0.2 GeV. Follow- 
ing Ref. PjHj], it is argued that right after these colli- 
sions, since the average pion separation is smaller than 
their interaction range, the pions in such a stage of the 
system evolution are still strongly interacting with each 
other. The effects of interaction among pions could then 
be modeled by considering that they are moving in an 
attractive mean field potential, which extends over the 
whole pion system. This implies, for instance, that in 
the two-pion case, they would not suffer any other ef- 
fects besides the mean field attraction and the identical 
particle symmetrization. Consequently, rather than be- 
ing in a gas, the pion system should be considered in a 
quasi-bound liquid phase, with the surface tension |l2j] 
acting as a reflecting boundary. Although details on this 
reflection depend on the pion wavelength, the pion wave 
function could be considered as vanishing outside this 
boundary. The pions become free when their average 
separation is larger than their interaction range. Due to 
the short range of the strong interaction, however, we 
would expect this liquid-gas transition to happen very 
rapidly, in such a way that the momentum distribution 
of pions could be essentially governed by their momen- 
tum distribution just before they freeze out. Under these 
circumstances, we would also expect that the observed 
pion momentum distribution would be modified by the 
presence of this boundary. And this is, in fact, what is 
analyzed in this work, as well as in Ref. R-|] 1 



On the other hand, since pion interferometry is sensi- 
tive to the geometrical size of the emission region as well 
as to the underlying dynamics, we would expect that the 
boundary would also affect the correlation function, but a 
priori we would not know how. Would it affect single- and 
double inclusive distributions similarly? How would be 
intercept of the two-particle correlation function behave? 
How would the general shape of this function be affected? 
For an insight into these questions, we here investigate 
the effects exerted by the boundary on the two-particle 



correlation function. We could naively expect that the 
importance of quantum statistics would progressively in- 
crease as the dimension of the emission region decreases. 
The results turned to exactly fulfill these expectations. 
Consequently, semi-classical approaches would have their 
applicability limited by the size of the emission region 
in focus. In other words, small emission volumes would 
stress the need for quantum statistics and, as a conse- 
quence, classical density matrices would lead to inconsis- 
tent results. This problem is clearly illustrated later in 
the present work. 

The plan of this paper is the following: in section II, 
we derive the single-inclusive distribution, as well as the 
two-pion correlation function, considering a density ma- 
trix suited for describing 7r ± 7r ± Bose-Einstein effects. In 
section III, the boundary effects on the two-pion corre- 
lation and single particle spectrum distribution are illus- 
trated by means of two specific examples. Section IV is 
devoted to illustrate the results which would be expected 
when previous methods for deriving two-pion interferom- 
etry formula are employed and the finite volume effects 
are studied. Finally, conclusions are discussed in section 
V. 



In momentum space, the corresponding pion creation 
operator, ^(p), and annihilation operator, VKp)> can be 
expressed as 



and 



where 



i£ + (p) = 5>j$(p) 



<Kp) = ^2a x ?p\(p) 



4>\(p) 



(27t) 3 / 2 



^A(x)e zpx dx. 



(5) 



(6) 



(7) 



We write the density matrix operator for our bosonic 
system as 



p = exp 



^{H - pN) 



where 



h = J2 e *«W> ^ = E«1 



A a A 



(8) 



(9) 



II. SPECTRUM AND TWO-PION 
CORRELATION FUNCTION 

In this section, we derive a generic formulation for 
the single- as well as for the two-particle inclusive dis- 
tributions, which would be suited for describing tt + tt + or 
7r~7r~ bounded in a finite volume. 

We assume the pion creation operator in coordinate 
space can be expressed as 



V?(x) = ]TaVl(x), 



(1) 



where a x is the creation operator for creating a pion in 
a quantum state characterized by a quantum number A. 
Then, V'a(x) is one of eigenfunctions belonging to a lo- 
calized complete set, which satisfies the orthonormality 
condition 

'dx^X^A'W^A.A', (2) 

and completeness relation 

5>A(x)V>A(y) = <5(x-y). (3) 

A 

Similarly, the pion annihilation operator in coordinate 
space can be written as 



V'(x) =^a A VA(x). 



(4) 



are the Hamiltonian and number operators, respectively; 
T is the temperature. 

The corresponding normalization is explicitly included 
in the definition of the expectation value of observables 
as, for instance, for an operator A 



{A) = 



tr{pA} 



(10) 



Then, the single-pion distribution can be written as 

Pi(p) = (^(pMp)) 

= EE^(p)v^(p)(«1«v) . (ii) 



A A' 



The expectation value (alay) is related to the occu- 
pation probability of the single-particle state A, N\, by 



(a\ay) = 6x,\>Nx 



(12) 



for a bosonic system in equilibrium at a temperature T 
and chemical potential p, it is represented by the Bose- 
Einstein distribution 



Nx = 



1 



exp [±(E X - p)] - 1 



(13) 



By inserting Eq. ( |12j ) and jl^ ) into (11), we obtain the 
single-particle spectrum for one pion species as 



Pi(p)=E JV ^a(p)^a(p). 



(14) 



The above formula coincides with the one employed in 
Ref. P] for expressing the single-pion distribution. 

Similarly, the two-pion distribution function can be 
written as 

P 2 (Pl,P2) = <^ t (Pl)^ t (P2)^(Pl)V'(P2)) 

= E V i A 1 (Pl)V^ 2 (P2)V'A3(Pl)0A 4 (P2) 

= E V'A 1 (Pl)'/ i A 2 (P2)V'A 3 (Pl)0A 4 (P2) 

Ai ,A2,A3,A4 



aA i aA 3 )(4 2 aA 4 )A 1 #A 2 



[' 

+ <«i 1 "A 4 )<aI 2 aA 3 )A 1 #A 2 

+ (« Ai "A 2 "A 3 OA 4 ) Ai = A 2 = A 3 = A. 

=Pi( P1 )i\( Pa )+x;E 



Ai A2 



^a 1 ^Va 2 V4(Pi)^a 1 (P2)^ 2 (Pi)V ; a 2 (P2) 
= Pi(pi)Pi( Pa ) + I ^7V a ^(Pi)^a(P2)| 2 . 



(15) 



Since we are considering the case of two indistinguish- 
able, identically charged pions, then 



(a^alaxax) = 2(a|a A ) : 



(16) 



From the particular form proposed for the density ma- 
trix in Eq. (8), we can see that (a A a A ) = (oa^a) = 0, 
showing that it would not be suited for describing 7r°7r° 
and 7r + 7r _ cases. For this purpose, the formalism pro- 
posed in Ref. f|,||,|ll| may be more adequate. 

The two-particle correlation can be written as 



C2(Pl,P2) = 



= 1 



Pi( P i)Pi(p 2 ) 

IEa^ a (Pi)^a(P2)| 2 



E a ^Va|V'a(pi)| 2 Ea^a|V'a(P2)| 2 



(17) 



We can see immediately from the above formula that 
if q = Pi — P2 = Owe have C2(p, p) = 2. We also 
notice that the result for the correlation function in Eq. 
( |l7| ) reflects the symmetrization over different states (and 
thus, the uncertainty in the determination of the pion 
state). 

Within this formulation we can also define the corre- 
sponding Wigner function, g(x, K), as 



^ K )=(2^E^/^( X +|)^( X -|) 



cxp(-iK • y)dy, 



(18) 



Consequently, we can write 
(^(pO^Pa)) = E^a(Pi)^a(P2) 

A 

= f e-ifa-P')-* 5 ( x ,k)dx. (19) 

By means of this Wigner function, the two-pion inter- 
ferometry formula can be re-written as |13-llq] 



C 2 (pi,Pa) = 1 



J e »q(x-y) g ( X;K)g(yiK)dxdy 
/ 5(x,Pi)g(y,p 2 )dxcZy 



(20) 



In the above equation, K = (pi + p 2 )/2 is the two- 
pion average momentum, and q = Pi — P2 is their rela- 
tive momentum. Here g>(x, K) can be interpreted as the 
probability of finding a pion at point x with momentum 
K. 



III. TWO-PION CORRELATION FROM A 
FINITE VOLUME 

A. Example 1 

In order to investigate the effect of the boundary on 
the single- and on the two-pion distribution functions, 
we assume that pions produced in high energy heavy-ion 
collisions are bounded in a sphere, just before freezing- 
out. In other words, their distribution functions are es- 
sentially the the ones they had while confined. The pion 
wave function should be determined by the solution of 
the Klein- Gordon equation 



[A + k 2 ] ip(r) = 



(21) 



where k 2 = E 2 — m 2 is the momentum of the pion. On 
writing the above equation, we have assumed that the 
potential felt by the pion inside the sphere is zero, while 
outside it is infinite. The boundary condition to be re- 
spected by the solution is 



V(r)| r =fl = 



(22) 



where R is the radius of the sphere at freeze-out time. 

The normalized wave function corresponding to the so- 
lution of the above equation can easily be written as 



V'fcZmW 



1 



RJ l+i (kR)\ r 




y, m (0,0)J,.i(*r) {r<R), 



(r > R). 

(23) 



The momentum of the bounded particle, k, can be de- 
termined as the solution of the equation 



J l+ i(kR) = 



(24) 



Inserting Eq. (23) into Eq. (7), we can determine the 
Fourier transform of the confined solution of a pion inside 
the sphere, as a function of the momentum p, as 



/2 1 / /? 

p* Ylm(p) RJ l+i (kR) [p^¥ 

P J l+i (pR)J lH (kR) - U l+ i( P R)J l+i (kR) 

(p ¥> k), 

i l Y lm (p) RJ l+ 3{kR) (p=k). 



(25) 



On deriving the above equation, we have made use the 
following integral equations 



I rdrJ l+1/2 (pr)J l+ i(kr) = 2 _ fc2 
P J l+ s(pR)J l+ ,(kR) - U l+k (pR)J l+i (kR)} , (26) 
and 

l rdrJ l+ r(kr)J lH (kr) = ^-Jf +i (kR) . (27) 



Besides, by imposing that the solution should vanish 
at the boundary, expressed by Eq. (|24|), it can be shown 
that 



lim 

p — >k 



pJi+3/2(pR)Ji+i/2(kR) - kJ l+ i(pR)J l+ i{kR) 



p2 _ fc2 



RJ*(kR) 



'i+i 



(28) 



i.e., ipklmip) is a continuous function of p at p = k. 

Using the condition expressed by Eq. (|I]), Eq. (|2E 
can be further simplified as 



J l+i (pR) . (29) 



i>ki m (p) = \j~i Yi m (p) 



p 



In terms of Eq. (|27]), the single- inclusive distribution 
function is given by 



p i(p) =^2N k i m ipi lm (p)i> k i m (p) 



kit, 



= E 



« m «p(^)-ivp 



l(-i) l Y? m (p)ffi 




-k- J l+ i{pR) 



p 2 — k 2 
2l + l\ fkJ l+ i(pRY 



(Eiipt) _ 1 \ 2np 



p 2 — k 2 



(30) 



where we have used that 

m— I 

X! Y hn(Pl) Y lm(p2) = 

m— — l 



21 + 1 

An 



Pi(pi-p 2 ). (31) 



Up to this point we have considered the pions confined 
in a sphere, which required the wave function to have a 
sharp change at r — R. However, as discussed in Ref. 
p6| , it could be more appropriate to consider a smoother 
boundary, by softening the potential felt by the pion at 
r = R. Unfortunately, this procedure would turn the 
problem into a very complex one |lq ], and beyond the 
scope of this paper. Nevertheless, as a diffuse boundary 
would cause a gradual decrease to zero of the pion wave 
function, it could be simulated by taking the limit R — > 
oo pHl in Eq. @, i.e., 



1 



Ji + i(pr)J l+ i(kr)rdr = -8{p-k) 



and 



J l+i (kR) 



~>Pklm = i L Yim(p) 



nkR 



pkR 



S(p-k). 



(32) 



(33) 



(34) 



Then, by imposing the completeness relation, Eq. (0), 
we can show that 



E^a(p)^a(p) = 



V 



(2tt)3 



(35) 



With the S function in Eq. (|3J) and Eq. (|3J), the 
single particle spectrum, in the limit R — > oo, is then 
written as 



Pi(p) 



1 



exp(^)-l 



V 



(2nf 



(36) 



where V — 4^i? 3 is the volume of the sphere. We see from 
Eq. ( p6[ ) that the ordinary Bose Einstein distribution is 
recovered in the limit of a very large volume. 
For p = and V = 4f R 3 , Eq. (BO) becomes 



Pi(p)| P =o = U^ 



^exp(^^)-l V 4 ^ 5 " 2 



where 



In the limit R — > oo, we have 



Pi(p) = 



V 

W) 



exp(^i) - 1 



(37) 



(38) 



(39) 



which is consistent with Eq. (p6|). 

From Eq. (B7|), we see that the intercept of the spec- 
trum depends on the value of the radius: as R increases, 
E n becomes smaller, and the maximum value of this dis- 
tribution, corresponding to |p| = 0, becomes higher. In 
all numerical estimates considered in the present work 
we have fixed /i = 0. In Figure 1, the normalized single- 
particle distribution is plotted as a function of |p|. We 
have chosen a discrete normalization, obtained by impos- 
ing Af — Yli=i Pi(j>i)t where N refers to the total number 
of bins in which the distribution is subdivided. 

We clearly see from Figure 1 that, due to the bound- 
ary effects, the maximum value of |p| in the spectrum 
decreases for decreasing volumes, being always smaller 
than the case corresponding to the R — > oo limit. The 
explanation for this behavior can be understood in terms 
of the uncertainty principle, i.e., as the volume of the 
system decreases, the uncertainty in the pion coordinate 
decreases accordingly, causing a large fluctuation in the 
pion momentum distribution. We should notice that this 
result coincides with the one obtained in Ref. |0,||, and 
is opposite to the results from Ref. ||0 . 

Regarding the spectrum, we could also inquire how 
would the freeze-out temperature affect it and how would 
the finite size effect compare with the R — > oo limit 
for different temperatures. This is illustrated in Figure 
2. The curves there should be compared in groups of 
two: solid (T = 0.14 GeV) and dotted (T = 0.11 GeV). 
For emphasizing the differences and similarities as |p| in- 
creases, we plot the difference between the two curves 
in each group, APi(p) = Pi(p)\n=3fm - i\(p)U=oo, in 
Figure 3. We see that the lower the temperature, the big- 
ger is the difference between the curves of each group, in 
the small momentum region. Decreasing the temperature 
has a similar effect on the spectrum as decreasing the ra- 
dius: in both cases the fluctuations in the small region of 
the pion spectrum increases and the corresponding max- 
imum is reduced. In other words, the boundary effects 
are more significant when we deal with systems whose 
dimensions and temperatures are small. 

We should observe that, except in Fig. 2 and 3 where 
the temperature dependence is studied, we have fixed 
T = 0.12 GeV. The reason for this relies on Shuryak's 
arguments fll2|| , according to which for temperatures in 
the range 0.1 < T < 0.2 GeV, the excited pionic matter 
would be better described as in a liquid phase inside a 
surface created by their mutual interaction. He added 
that, for T > 0.15 GeV, the influence of resonances be- 
come important but they are not included in the present 
study. Therefore we chose T = 0.12 GeV, which is also of 
the order of the recent experimental freeze-out tempera- 
ture estimated from both interferometry and spectra. 

Similarly, we can write for the expectation value of 
the product of two pion creation operators in momentum 
space 



^klmiVX^klmiVi) 



(^(p 1 )^(p 2 )) = ^. , 

H^ exp(^i) I 



£ 



i 



kim exp 

"2 

Pi 



(ta) 



H)'^(p0 



(iYY*m(p2) 



vi - 
-k 



P2 



J l+ i(piR) 

J 1+ i(P2R) 



£ 



ki exp 



(^) 



1 V PlP2 



Ji+dpiR) 



( P j-k3)(pi-kY l+ 



(40) 



The two-pion interferometry correlation function can 
then be estimated by inserting the above expression into 
Eq. (|l7|). In general, this function will depend on the 
angle between pi and P2- For the sake of simplic- 
ity, we will consider pi parallel to P2, implying that 
P/(pi • p>2 = ±1) = (±1) ■ The results for two-pion in- 
terferometry corresponding to different values of the pair 
average momentum K = (pi + P2)/2 but fixed temper- 
ature are shown in Figure 4. We can see that, as the 
pair average momentum increases, the apparent source 
radius becomes bigger, which is an interesting behavior, 
if we compare to results corresponding to expanding sys- 
tems. In this last case, the probed part of the system 
decreases with increasing average momentum. Naturally, 
our present approach does not consider the effects of ex- 
pansion and the enlargement of the system's apparent 
dimensions with increasing K, seen in Figure 4, has its 
origin in the strong sensitivity to the dynamical matrix. 
This can be better understood by observing the presence 
of the weight factor N\ in Eq.([17|), with N\ expressed 
in Eq.(|13|). The increase of the average momentum re- 
flects the increase in the individual momenta pi and P2, 
which comes from larger values of the sum coefficient, 
k, in Eq.dri]). This has two opposite effects: the factors 
l/(p| — k ) give larger contribution for k ~ pi. However, 
bigger values of k would also make the exponential factor 
(with [i = 0) drop faster. Being so, by increasing K we 
are effectively including a larger number of k coefficients 
that contributes to the sum in Eq.(p7|), with decreasing 
weight ~ exp (—Eki/T). The interference of these extra 
terms corresponding to larger k with the terms already 
considered in the sum would make the correlation func- 
tion drop faster, consequently becoming narrower. Al- 
ternatively, we could understand these results by notic- 
ing that pions with larger momentum come from larger 
quantum A states which, in turn, correspond to a smaller 
spread in coordinate space. As the weight factor in Eq. 



( |17| ) is of Bose-Einstein form, larger quantum states will 
give a smaller contribution to the source distribution, 
causing the effective source radius to appear larger. In 
order to confirm that the weight factor in Eq.(|17|) is the 
responsible for the behavior observed in Fig. 4, let us 
consider the case in which we choose it to be a constant 
factor, for instance N\ — 1. This situation could be de- 
rived from the Bose-Einstein distribution form by consid- 
ering T 3> 1, so that the two-pion interferometry results 
would become insensitive to the average momentum, due 
to the very large values of the temperature. The numeri- 
cal result corresponding to this case is also shown in Fig. 
4 (narrower curve). On the other hand, with the help 
of the completeness relation, Eq.(||), and of Eq.(|7]), by 
also assuming that the pions are confined in a sphere, it 
is straightforward to derive the following K independent 
two-pion correlation function 



C 2 (q) 



1 



9 , sin(o_R), 

'Rcos(qR) ^-^] 2 . (41) 



q 4 R 6 



For completeness, we also include in Fig. 4 the curve 
based on Eq. (|4l|) , which coincides with our numerically 
generated curve, cross checking the correctness of our 
numerical calculation. 

Figure 5 shows the two-pion correlation function for 
increasing values of the spherical radius, i.e., for enlarging 
volumes. From it, we can clearly see that, as the confining 
volume increases, the source radius derived from two-pion 
correlation also increases, as would be expected. 

Again, as discussed previously for the spectrum, we 
could estimate the effect of a diffuse boundary on the two- 
pion correlation function by considering the limit R — > 
oo. By inserting Eq. (134) into Eq. (Boh, remembering 
that in this limit we can take ^2 k — ► J dk, and using the 
previous result for the spectrum in this limit, Eq. (36), 
we finally obtain that 



C 2 (pi,p 2 ) 



(Pi ^ P 2 ) 
(Pi = P2) 



(42) 



as would be expected. 

To conclude this section we should keep in mind that, 
if the system size is very small, it would be sensitive 
to the boundary effects even if we considered a diffuse 
boundary. On the contrary, if the system size is very 
large, we would not expect a significant effect in neither 
the sharp nor the diffuse boundary case flq] , 

B. Example 2 

In this subsection, we will study the sensitivity of spec- 
trum and of the two-pion correlation function to the sys- 
tem boundaries, by considering the pion system inside a 
box of dimensions L x L x L. We choose first periodic 



boundary conditions. In this case, the eigenfunction can 
be written as 



^k(r) = -=exp(-ik-r) 



(43) 



Here k is the quantum number which satisfies following 
constraint 



ki ■ L = 2 Hi TT 



1=1,2,3 

n; = 0,±l,±2,. 



(44) 



The corresponding Fourier transform V'k(p) can be ex- 
pressed as 



"0k (P) = 



1 



N 



(2^)3/2 yy 



sin[(pi-Ai)L/2] 



Pi - ki 



sin[(p 2 -fc2)L/2] 



2'2 



sin [(p 3 - h)L/2] 



P3 



(45) 



The single-particle and two-particle distributions fol- 
low from Eq. (14) and (|15|). We should notice that, in 
the limit L — > 0, and using the condition (H), we find 
the contribution of only one state (k = 0) to the two-pion 
correlation function, resulting in 



C(pi,p 2 ) = 2 



(46) 



On the other hand, if we take the limit of V — > 00, Eq. 
tBI) becomes 



^k(p) 



1 



(27r) 3/2 (5(p-k). 



(47) 



With the above form for V'k(p) in the limit of very 
large volumes, we obtain for the correlation function 



C 2 (P1,P2) 



(Pi + P2) 
(Pi = P 2 ) 



(48) 



If, instead of the periodic boundary conditions, we con- 
sider that the pions are confined in the box, i.e., we as- 
sume the potential outside it is infinite, then two classes 
of solutions are possible: 




sin(fci ■ x) sin(fc2 • y) sin(fc3 • z) 



1=1,2,3 

m = 1,2,. 



(49) 



(50) 



and 

^ /7 (x) = \j yj cos(fci • x) cos(k 2 ■ y) cos(fc 3 • z) , (51) 
with 



fcj • L = (2n, — 1)tt 



1=1,2,3 
n-i = 1,2,. 



(52) 



It can be shown that, for V — > 0, we have 
C 2 /J/ (Pi,p 2 ) = 2, 



while, for V — ► oo, we obtain 



C2' 7/ (pi,p 2 ) = | 2 



(Pi + ±p 2 ) 
(Pi = ±Pa) 



(53) 



(54) 



The reason for including the (±) signs in Eq .(p4j) comes 
from the parity property of Eq. ( |49|) and (glj). From 
them, it is immediate to see that solution I has negative 
parity, while solution II has positive parity. The corre- 
sponding Fourier transforms then show the same parity 
property, i.e., 



&(p) = -&(-p) ; & z (-p) = & j (p) 



(55) 



From the above results, we can show that the single 
particle spectrum and two-pion correlation function cor- 
respondingly have the following properties 



rjl.lll \ r>I,II/ \ 

A (p) = -Pi (-p), 



iiji 



i, ii t 



C 2 ' (Pi,-Pa) =C 2 ' (Pi.Pa) 



(56) 



(57) 



and 



In particular, we see from Eq. (p3J), (|54|), and (p7 
that, if we choose Pi = p 2 = p, we immediately get 
C 2 ' (pi = p 2 = p) = 2 for the confined boundary con- 
dition in both volume limits. That is the reason why, as a 
consequence of the parity property of the wave function, 
Eq. ( |54J) could be extended to C 2 ' (pi = ±p 2 = p) = 2. 

For periodical boundary condition, however, we have 



^k(p) = ^-k(-p) 



(58) 



Then, for the single particle distribution we will have the 
following relation 



Pi(p) = Pi(-p) 



(59) 



Nevertheless, the two-pion correlation function, which 
can be written as 



C 2 (Pi,-p 2 ) = 1 



IEk^k( P i)^k(- P2 )| 2 
Ek^k|4(pi)| 2 E k «k(p 2 )| 2 ' 

(60) 



in the case of periodical boundary condition, will show 
no well-defined property under momentum reflection as 
the one expressed by Eq. (p7|). 



IV. CONVENTIONAL HBT FORMULATION 

We now discuss the case of the conventional formula- 
tion, usually adopted in HBT analysis, in terms of clas- 
sical currents j(x) p] representing the pion sources. For 
simplicity, we consider the momentum as the only quan- 
tum number involved in the problem, i.e., we denote {A} 
as {p}. Besides, we also assume that the pion state could 
be characterized by the measured momentum. For in- 
stance, V'(Pi) represents a pion in a quantum state de- 
noted by pi . Then the single-particle and the two parti- 
cle distributions can be written as 



Pi(p) = jvw(pMp), 



(61) 



and 



P 2 ( P i,p 2 ) = 7V( P i)iV(p 2 )^(pi)r(p 2 )V(Pi)V'(p 2 ). 

(62) 

In the above relations we dropped the subscript {A} 
of the state. N(p) is the Bose-Einstein distribution. In 
order to connect to the HBT effect, we need to make one 
further assumption: we assume that the source j(x) is 
chaotic and a function of the coordinates only. ip(x) is 
determined as a solution of the equation 



(A+p 2 W(x)=j(x) 



(63) 



The superscript j is introduced as a reminder that if) 3 
is the solution of equation (|63|), in the presence of the 
source j(x). Then ij) 3 {x) can be written by 



i> j (x) = / G(x,x')j{x')dx' 



(2tt) 3 / 2 



-ip-(x-x') ■( i 



i(x')dx'. 



(64) 



In the above expression we have used the fact that 
j(x') is localized in a small volume. The corresponding 
function in momentum space is then 



(65) 



^ J (P)= / j(x)e lpx dxi. 
The current j(x) can be expressed as 

N 



(66) 



i=l 



where i denotes the number of the collision center; Ai 
is a weight factor which represents the amplitude of the 
emitter. Assuming they are chaotic and a function of the 
coordinates only, then 



{j*(x)jj(y)} = $ij3i{x)ji{y)- 



(67) 



Here, {■ • •} denotes average over phases. The single 
particle spectrum and two-pion distribution function are 
then written as: 



P 1 (p)=N{p)J2\M%(P)\ i 



(68) 



and 

P 2 (Pi,P2) = N( Pl )N(p 2 ) 



k 



\A k \ 2 \j k (Vit 



N 

IE 

i=l 

N 

~E 



\A\ 2 j*(Px)ji(P2)f 



\A i \ A j:(pi)jUP2)j l (Pi)MP2)}, (69) 



with 



3i(p) 



(2tt) 3 / 2 



3i{x)^dx. 



(70) 



Inserting Eq. (|68| ) and (|69| ) into (|17|), we obtain the 
two pion intcrfcrometry formula expressed as 



C 2 (pi,p 2 )= 1 



EJ^| 2 J*(Pl)ji(P2)| 5 



E l l^l 2 |j I (pi)l 2 EJ^I 2 b fc (P2)l 2 
EJ^I 4 bT(pi)^(P2)l 2 



EJ^I 2 lj,(pi)l 2 E fc l^l 2 lj & (P2)l i 



(71) 



According to Ref. [g| ,the strength of each current could 
be localized around some inelastic scattering center Xi, 
such that 



3i(x) = j(x-Xi) ; \A,\ 2 = p(xi) 



(72) 



The current j(x— xi) is considered to be peaked around 
Xi, and could be characterized by the size scale of the 
wave packet; p(xi) is the source distribution function of 
the emitter. Naturally, we are now considering a simpli- 
fied picture, in which phase-space correlations are absent. 
Then Eq. ( J7l| ) can be further simplified as 

Ei,j p{x l )p{x j )cos[{pi -p 2 )(xi - Xj)] 



C 2 (pi,P2) = 1 



Ei P( x i) Efc P( x k) 
T,iP(Xi)p(Xi) 



T, t p(^)T, k p( x k)' 



(73) 



The last term in Eq. @, (|l|) and @ discounts the 
contribution corresponding to emitting two pions from 
the same source. In the case of very large volumes this 
type of contribution is usually considered to be small 
when compared to the emission from separate sources. 
Consequently, in cases where this term O (1/V) can be 
neglected [p 15 , we recover the well-known two-pion in- 
terferometry formula 



C2(Pi,P2) = l + p{x)p(y)cos[(p 1 ~p 2 )(x~y)]dxdy. 

(74) 

In Ref. JL4J , a general semi-classical approach to two- 
pion interferometry was used, in which a Gaussian wave 
packet spread was allowed for incorporating minimal ef- 
fects due to the uncertainty principle. As pointed out in 
the above reference, Eq. (f74|) corresponds to approaching 
the classical regime, i.e., it would be valid only when the 
wave packet size is negligible, which would also be equiv- 
alent to consider system sizes much bigger than the wave 
packet size. Being so, the above derivation would be con- 
sidered as a good approximation only in cases where the 
source size is large, as in heavy-ion collisions. However, 
we should be cautious when using it in e + e~ collisions 
as, in that case, the source radius is small and the third 
term in Eq. ([73]) may not be negligible. Besides, the 
chaotic source ansatz is also questionable there. Just to 
emphasize this point, let us naively consider a fictitious 
source of fm size, i.e., p{x) — 5(x). Then, from Eq. 
(f74|), we would get 



C 2 (pi,p a ) = 2, 
and, for the chaoticity parameter 

A = C 2 (p,p)-l = l. 



(75) 



(76) 



Naturally, we cannot think of a "zero size" source as 
being chaotic. Actually, as it has been stated in Ref. [p), 
and illustrated above, the chaotic ansatz is only correct 
when the source size V is much larger than the size of 
the wave packet. In the fictitious source case above, if 
we do not neglect the third term in Eq. (|7l|), we would 
obtain 



C 2 (Pi,P2) = 1, 



(77) 



which is the correct two-pion interferometry result for the 
above model p7| . 

For confronting the role of the correction term in Eq. 
( |7l| ) and Eq. (73) with the correlation function estimated 
by using Eq. (|74| ) we choose, for simplicity, A; = 1/VV 
and ji(x) — 1/^/Va, where V is the total volume of the 
source and V a is the volume of the emitter, which is of the 
order of the wave packet size. In Figure 6 we show the 
corresponding results. We see that, for V — (5/m) 3 and 
V a = (1/to) 3 , the third term is much smaller then the 
second, and could be neglected, while for V — (2/rn) 3 
and V a = (1/to) 3 , its correction is substantial. Clearly 
then, these results depend strongly on the wave packet 
size: the smaller it is with respect to the system size, 
the better is the approximation represented by Eq. (174) . 
From the results in Figure 6, it seems that, if the wave 
packet size is about 1/to, we could use the conventional 
pion interferometry formula in Eq. (f74[) for analyzing 



pion interferometry in heavy-ion collisions. However, 
we could not use it to analyze e + e~ collisions, since in 
that case the source radius is of the same order as the 
wave packet size, and the contribution the third term is 
non-negligible. It is interesting to notice that the above 
derivation is equivalent to the one in Ref. Iqjlq] where a 
density matrix formulation was also used. Appendices I 
and II contain, respectively, further discussion regarding 
a density matrix formulation leading to an equivalent of 
Eq. (|7l|), and a simple unified form for the two formula- 
tions. 



V. CONCLUSIONS 

In this paper, we derive the two-pion correlation func- 
tion by adopting a different density matrix, as given in 
Eq.(||). The finite volume effects on the pion spectrum 
were then studied in Figures 1, 2 and 3, leading to sim- 
ilar results as in Ref. |7p|JlO||. We found that the small 
momentum region is depleted with respect to the Bose- 
Einstein distribution. The effects on two equally charged 
pion correlation function were also analyzed. The results 
in Figures 4 and 5 show that the correlation function 
shrinks for increasing average pair momentum and for 
increasing size of the emitting source, respectively, corre- 
sponding to an increase of its inverse width. The first re- 
sult reflects a strong sensitivity to the dynamical matrix, 
through the Bose-Einstein weight factor, as discussed in 
Section III. A. We also discussed the effects of a diffuse 
boundary on the spectrum and correlation function by 
considering a smooth decrease to zero of the wave func- 
tion as R goes to infinity. 

We compared as well the results obtained by means of 

' or dH, 



Eq. (|17| ) with those estimated by using Eq. ( |71|) 
showing that they may differ significantly when small vol- 
umes are considered. For instance, from Figures 4 and 5, 
and as a result of Eq. ( |l7|) , we see that the boundary af- 
fects the single- and the double-inclusive distributions in 
a consistent way, as the intercept of the correlation func- 
tion remain unchanged when altering the system size. 
Nevertheless, when looking at the curves generated by 
using Eq. ( f71[) in Figure 6, we see that the intercept of 
C{q) drops as the size of the system is reduced, leading to 
unphysical results. The relation of the system size to the 
wave packet size is the key to understand this behavior: 
when this last one cannot be considered much smaller 
than the first one, the additional term in Eq.(j7l|) or 
Eq.([73|) would give a non-negligible contribution. Other- 
wise, Eq.(|7l|) or Eq.(j73]) would approximate the conven- 
tional HBT formulation, represented by Eq.(|74|), which 
was derived under the condition that the system size is 
much larger than the wave packet size. Nevertheless, 
we should notice that this assumption was not necessary 
for obtaining Eq. flT7|), since this was derived strictly 
within the quantal realm. It is interesting to remark 



that, in the simplistic case of a Gaussian breakup dis- 
tribution, considering wave packets with non-negligible 
widths, Eq. (61) of Ref. [[Uj showed that the inverse 
width is enlarged in direct proportion to the wave packet 



size, R\ = R 2 



Aa 



In this expression, R 



4(P 2 +Ap 2 )- 

is the Gaussian width in space-time, P is the width of the 
Gaussian distribution in the momentum space; Ax and 
Ap are the corresponding wave packet spread. For R ~ 2 
fm, P ~ 0.14 GeV/c, and minimum packets with Ax = 1 
fm, the corresponding inverse width would be R/\ ~ 2.16 
fm, an increase of about eight percent. This very rough 
estimate seems to be of the same order of the decrease 
in width (i.e., increase in the apparent radius) seen in 
Figure 6. Finally, in Appendix I we show that, by adopt- 
ing another density matrix instead of the one proposed 
in Eq. (B), we can derive an interferometry result which 
is similar to Eq. ([71J). And in Appendix II, we suggest a 
simplified way of unifying these two formulations, i.e., of 
recovering each of them by means of a parameter choice. 
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VI. APPENDIX I 



VII. APPENDIX II 



In what follows, we show that, by considering the fol- 
lowing density matrix 



P 



n* 



(78) 



instead of the one proposed in Eq. (||) , we can also derive 
an interferometry result which is similar to Eq. (71), 
where 



Pa = 



E 



(ol) B * ... exp(-n A • (E x - Hx)/T) (a A ) r 



'n\\ 



-|0) 



n\\ 



'n\\ 



(79) 



The corresponding pion multiplicity distribution is 
given by 



Pn 



{nx 



n x < 



exp(-(n A )), 



(n x ) = N x = exp(-(E x -n)/T), 
and the single particle spectrum is then 

P!(p)=^iV A |^(p)| 2 . 



(80) 



In the limit V — ► oo, we obtain the Boltzmann distri- 
bution instead of Bose- Einstein distribution derived pre- 
viously in Section III. 

Similarly, we can derive the two-pion correlation func- 
tion as 



C 2 (pi,P2) = 1 



\ExN^l(pi)MP2)\ 2 
ExNx\Mpi)\ 2 ExN x \Mp2)\ 2 

J2 x n 2 x \Mpi)Mp2)\ 2 
E x n x \Mpi)\ 2 J2 x n x \Mp2)\ 2 ' 



(81) 



The major difference between this equation and Eq. 
(O) has its origin in the fact that now we have 



{a\a x a x a x ) 



(4«a) 2 



(82) 



instead of Eq. (pit), in section II. This comes from the 
fact that, in Eq. (|82j), we are dealing with classical (i.e., 
distinguishable) particles. In Appendix II, we show a 
way to unify the formulation leading to Eq. ( |17[ ) and Eq. 

©■ 



If we assume that a single pion state could be described 
by the quantum number A, then the single particle spec- 
trum could be expressed as 



Pi(p) = 5>a|&(p)| s 



(83) 



Here cu x is the occupation probability of the single- 
particle state A. In the two particle distribution case, 
the pions could be described by two different quantum 
numbers, Ai and A2. By imposing the symmetrization re- 
quired by the Bose-Einstein statistics, the two-pion wave 
function could be written as 



1 



V%,A 2 (Pl,P2) = — =[V'A 1 (Pl)'0A 2 (P2) + 
V'A 1 (P2)^A 2 (Pl)]- 



(84) 



If the two-pions are in the same state, then we would 
have 



^a,a(Pi, P2) = A ■ ^ a (pi)-0a(P2)- 



(85) 



If the multiplicity distribution follows the geometry 
distribution, then A = y/2. If, however, multiplicity dis- 
tribution has Poisson form then A = 1 . Being so, the 
two-pion distribution could be expressed in terms of A 
as 

-P2(Pl,P2)= ^2 W A!^A 2 

Ai,A2,Ai^A2 

i[bk( P l)| 2 |^A 2 (p2)| 2 + 

hM Pl )| 2 |<MP2)| 2 + 

V^ Al (Pl)^A! (P2)^A 3 (P2)V% (Pi) 
V^ Al (Pa)^! (Pl)^A2 (Pi)^a 2 (Pa) 

+\A\ 2 Y,u 2 x \Mpi)i>x(P2)\ 2 

A 
= Pl( P l)Pl(p 2 ) + 

|^W A ^( P l)^ A (p 2 )| 2 
A 

-2^^|^ A ( Pl )^ A (p 2 )| 2 

A 

+ |A| 2 ^ w 2 |^ A ( Pl )^ A (p 2 )| 2 (86) 



Consequently, for different choices of A, i.e., A = \/2 
or A = 1 discussed above, we could recover the results in 
Eq. (H|) or Eq. (|l|), respectively. 
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FIG. 1. Normalized spectrum (in arbitrary units) vs. mo- 
mentum |p| (in GeV/c). The input temperature is T = 0.12 
GeV and the chemical potential is [s, — 0. The solid line cor- 
responds to the Bose-Einstein distribution, i.e., to the limit 
R — * co. The dotted and dashed lines correspond, respectively, 
to the R = 6 fm and R = 3 fm cases. 
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FIG. 2. Normalized single-particle distribution (in arbi- 
trary units) vs. momentum |p|. The input parameters are 
(i = 0, and the radius in the finite case, R = 3 fm. The 
curves are shown for two values of the freeze- out temperature: 
T = 0.14 GeV (solid) and T = 0.11 GeV (dotted) and com- 
with the corresponding ones in the R — > oo limit. 
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FIG. 3. Difference in the momentum distribution (in arbi- 
trary Units), APl(p) = Pl(p)|fl=3/m - -Pl(p)|fl=oo, vs. |p|, 

of curves with R = 3 fm and the corresponding ones in the 
R — > oo limit, for two values of the freeze-out temperature T, 
as indicated in the plot, with fi = 0. 
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K = 0.3 GeV/c 
- Eq.(41) 

numerical results 
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FIG. 4. Two-pion correlation function, C^q), is shown 
versus the momentum difference |q|=|pi — pa|. The input pa- 
rameters are T = 0.12 GeV, R — 3 fm, and [i — 0. The solid, 
dotted, and dashed lines correspond to the average pair mo- 
mentum values K— 0.1 GeV/c, 0.2 GeV/c, and 0.3 GeV/c, 
respectively. The circles refer to numerical results similar to 
the previous ones but with unity weight factor, N\ — 1, in Eq. 
nil). The dot-dashed line corresponds to the analytical result 
in Eq. (41). 
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FIG. 5. The correlation function, 62(5), is plotted as a 
function of the momentum difference, |q|=|pi — p?\. The in- 
put parameters are the temperature T — 0.12 GeV, the average 
pair momentum K — 0.4 GeV/c, and /j, — 0. The solid, dot- 
ted, and dashed lines correspond, respectively, to the sphere 
radius R = 2 fm, 3 fm, and 7 fm. 
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FIG 



6. Two-pion correlation function, C2(q), is shown 
|q|— |Pi — Pal- For helping visualization the curves are 
by brackets in two groups. In the first, signaled by 
(k), the curves were obtained with the help of Eq. (|7ij). The 
other group, signaled by (#), corresponds to curves obtained 
by means of Eq. n?A). The values adopted for the finite sys- 
tem sizes, i.e., L= 2 fm, 3 fm, and 5 fm, are shown in the 
plot. 
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